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Control  of  an  Unstable,  Nonminimum  Phase  Hypersonic  Vehicle 


Model 

Michael  W.  Oppenheimer  * 

David  B.  Doman  t 

Air  Force  Research  Laboratory,  WPAFB,  OH  45433-7531 

Abstract 

In  this  work,  a  control  law  for  an  unstable,  nonminimum  phase  model  of  a  hypersonic 
vehicle  is  developed.  The  control  problem  is  difficult  due  to  the  locations  of  the  plant  poles 
and  zeros.  For  an  unstable  system,  feedback  is  required  to  stabilize  the  plant.  However,  one 
cannot  make  the  loop  gains  arbitrarily  large  without  driving  one  or  more  of  the  closed-loop 
poles  into  the  right-half  of  the  s-plane,  since  the  system  is  nonminimum  phase.  Thus,  there 
is  a  limited  range  of  feedback  gain  that  results  in  a  stable  system.  The  nonminimum  phase 
behavior  also  places  restrictions  on  the  closed-loop  bandwidth.  For  the  hypersonic  vehicle 
control  problem,  low  frequency  control  is  desired  and  a  rule  of  thumb  is  that  the  closed- 
loop  bandwidth  must  be  less  than  one-half  the  right-half  plane  zero  location.  A  right-half 
plane  zero  located  in  the  region  of  the  desired  gain- crossover  frequency  makes  it  impossible 
to  achieve  the  desired  level  of  tracking  performance.  The  achievable  closed-loop  bandwidth 
might  be  so  small  that  adequate  control  of  the  system  is  not  achieved.  Direct  cancellation 
of  the  right-half  plane  zero  with  an  unstable  pole  in  the  controller  is  not  an  option.  In  this 
work,  a  modified  dynamic  inversion  controller  is  developed  for  a  linear,  time-invariant  model 
of  a  hypersonic  vehicle.  This  modified  dynamic  inversion  controller  differs  from  the  standard 
dynamic  inversion  approach  in  that  it  does  not  attempt  to  cancel  the  right-half  plane  zero 
with  a  pole,  instead,  it  retains  right-half  plane  zeros  in  the  closed-loop  transfer  functions  and 
uses  an  additional  feedback  loop  to  stabilize  the  zero  dynamics. 

HSV  Model 

The  model  to  be  controlled  in  this  case  is  given  by 

x  =  Ax+Bu 
y=Cx+Du 
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where  A  G  Mn:m,  B  G  Rnxm ,  C  G  D  G  Rpxm,  and  n  =  9,  m  =  p  =  2.  The 

model  outputs  are  velocity  (ft/sec)  and  flight  path  angle  (deg),  while  the  inputs  are  elevator 
deflection  (deg)  and  temperature  addition  in  the  combustor  (deg  R).  More  specifically,  the 
state-space  matrices,  derived  from  a  nonlinear  model,1  are 


A  = 


— 4.8e-4 

2.05 

0 

— 5.1e-6 

-32.17 

.084 
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0 

— 5.8e-7 
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0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

1 

0 

0 

-.068 

-7368.3 

0 

.0159 

0 

-286.5 

-.641 

1.518 

1.078e_: 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-.076 

-5668.22 

0 

.018 

0 

-4.72 

8.82e-4 

-401.39 

-.78 

(2) 


B  = 


-62.589 

-.0222 

-0.997 

0 

0 

0 

139.032 

0 

-2500.6 


0.0064 
— 2.61e-7 
8.25e-7 
0 
0 
0 

-.000115 

0 

-.0769 


(3) 


C  = 


1  0  0000000 

0  -1  0010000 


(4) 


where 


D 


0  0 
0  0 


(5) 


x 


V  OL  Q  h  0  T)f  f)f  T)a  f)a 


1  T 


(6) 


V  denotes  velocity  or  true  airspeed,  a  denotes  the  angle  of  attack,  q  denotes  the  pitch  rate, 
h  is  the  altitude,  0  is  the  pitch  attitude,  and  the  last  four  states  represent  temporal  modal 
coordinates  that  describe  the  first  bending  mode  of  the  fore  and  aft  fuselage.  The  outputs 
of  the  model  that  we  wish  to  control  are  the  velocity  V  and  the  flight  path  angle  7  =  6  —  a. 
The  idea  in  this  work  is  to  use  a  dynamic  inversion  type  scheme  to  control  the  plant  and 


2  of  18 


provide  a  desired  response.  Unfortunately,  this  plant  is  nonminimum  phase,  with  poles  and 
transmission  zeros  at  the  following  locations: 


Poles  =  {-0.39  ±  j'20.03,  -0.32  ±  jl6.94,  -1.0035,  0.9322,  -5.6e-4,  -3.88  ±  j'0.041} 
Zeros  =  {0, 1.949,  -1.948,  -0.391  ±  jl9.58,  -0.321  ±  jl6.95} 


(7) 


Hence,  dynamic  inversion2  (DI)  is  not  an  option  due  to  the  right-half  plane  zero.  This  is 
because  DI  will  try  to  cancel  plant  zeros,  thus  resulting  in  a  right-half  plane  pole  in  the 
closed-loop  system.  Fortunately,  a  DI  type  scheme  can  be  used  to  circumvent  this  issue. 
Dynamic  extension3  retains  the  nonminimum  phase  zero  (does  not  try  to  cancel  it  with 
an  unstable  pole),  yet  still  has  the  desirable  trait  of  decoupling  the  system.  The  following 
describes  the  development  of  a  dynamic  extension  algorithm  for  the  system  described  above. 


Dynamic  Extension 

Consider  the  system  described  in  Eq.  1.  The  vector  y  is  a  set  of  controlled  variables 
(CVs).  Dynamic  inversion  is  used  to  decouple  the  system  and  produce  desired  responses 
from  the  CVs.  To  begin  the  development,  consider  each  CV: 

Hi  =  C,x  f  D,u.  i=  1.2  (8) 


Since  each  Dj,  i  —  1,2  is  a  row  of  zeros,  the  input  does  not  explicitly  show  in  the  output 
equation.  To  obtain  an  equation  with  the  input,  differentiate  y  to  obtain 

Vi  =  C  ii  =  Ci  Ax  +  QBu  (9) 


In  this  case,  QB  ^0,*  =  1,2,  so  there  is  no  need  to  take  any  more  derivatives  of  the  output. 
Now,  the  system  dynamics  are  given  by 


x  =  Ax+Bu 
y  =  CAx+CBu 


(10) 


In  order  to  decouple  the  controlled  variables,  the  matrix  CB  must  be  right  invertible.  If  this 
is  the  case,  then 


u  =  (CB)-1  (ydes  -  CAx) 


(11) 
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where  ydes  is  a  vector  of  pseudo-controls.  In  this  case, 


Yd 


es 


Vtdes 

^des 


(12) 


where  Vt  is  the  vehicle’s  velocity  and  7  is  the  flight  path  angle.  Applying  the  vector  u  in 
Eq.  11  to  the  output  and  state  dynamics  equations  produces 


y  =  CAx+CBu  =  CAx  +  (CB)  (CB’1)  (y,„  -  CAx)  =  yde,  (13) 


so  that  the  controlled  variables  follow  exactly  the  pseudo-controls  and 

x  =  Ax+Bu  =  Ax  +  B  (CB-1)  ( ydes  —  CAx) 

=  (A-B(CB)-1CA)x  +  B(CB-1)y,„  (14) 

—  A*  +  B  (CB-1)  yto 

As  it  turns  out,  the  eigenvalues  of  Ax  are  the  poles  of  the  zero  dynamics  and  are  identically 
equal  to  the  zeros  of  the  original  system  (see  Eq.  7).  So,  the  state  dynamics  are  determined 
by  the  inverse  state  space  system 


x  =  A'1'  B  (CB  ')  v 
x=Ix+0v 


(15) 


Let  S  be  a  matrix  of  left  eigenvectors  of  Ax.  Then,  Eq.  15  can  be  transformed  into  the 
following  Jordan  form 

£  =  A£  +  SB  (CB-1)  v  (16) 

where  A  =  SA^S-1.  The  states  in  Eq.  16  are  related  to  the  original  states  through  the 
transformation  £  =  Sx. 

The  transmission  zeros  are  defined  as  the  values  of  £,  the  vectors  zj,wi ,  the  input  zero 
directions  or  zo,wq,  the  output  zero  directions  that  satisfy 


- 1 

hH 

1 

> 

l 

W 

_ i 

z  I 

[  CA  CB 

W  / 

L 

J  L 

j 

1“ 

- 

r\ 

l—H 

1 

> 

-B 

z  o  Wo 

CA 

CB 

(17) 


(18) 
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From  the  output  zero  case,  solve  for  z0  to  get 


— z0B  +  w0  (CB  x)  =  0 


(19) 


Solving  for  w0  yields 


w0  =  z0B  (CB  x)  (20) 

Substituting  Eq.  20  into  the  first  output  equation  in  Eq.  18  yields 

z0  (Cl  -  A)  +  w0CA  =  z0  (Cl  -  A  +  B  (CB"1)  CA)  =  0  (21) 

Hence,  the  transmission  zeros  are  the  eigenvalues  of  the  matrix  (A  —  B  (CB-1)  CA),  which 
is  identically  equal  to  the  matrix  Ax  defined  above.  This  result  shows  that  the  transmission 
zeros  of  the  original  system  are  poles  of  the  zero  dynamics  of  the  dynamic  inversion.  In  order 
to  obtain  stable  zero  dynamics,  the  desired  dynamics  of  the  CVs  must  be  modified.  This 
requirement  is  due  to  the  nonminimum  phase  behavior  of  the  original  system.  Additionally, 
the  vectors  z0  are  the  eigenvectors  of  Ax .  Gathering  all  these  eigenvectors  into  a  matrix,  Z0, 
it  is  seen  that  Zo  is  a  matrix  of  left  eigenvectors  of  Ax.  Hence,  Zo  =  S.  Then,  from  Eq.  20, 

W0  =  Z0B  (CB-1)  =  SB  (CB-1)  (22) 

Modified  Dynamic  Inversion  Controller 

Case  1 

Suppose  that  one  transmission  zero,  <^,  is  in  the  right-half  plane  and  would  thus  result  in 
unstable  zero  dynamics.  Without  loss  of  generality,  suppose  that  the  output  zero  direction 
is  given  by 


w 


RHP 

0 


0  w. 


RHP 


(23) 


This  means  that  only  pseudo-input  2  excites  the  bad  zero  state  (the  state  corresponding  to 
the  nonminimum  phase  zero).  Hence,  we  select  i)des2  to  include  a  stabilizing  term,  from  the 
bad  zero  state  (£,RHP),  such  that 


=  ides  +  kiiRHF 


(24) 
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In  Jordan  form,  the  dynamics  of  the  bad  zero  state  become 


cRHP  .  „trhp 


—  Ci£RHP  +  OUi  +  W0H2V2  —  C i£RHP  +  w0il2  (' jdes  +  k^RHP^(25) 


Thus, 

cRHP  wQHP2  Wes 

^  8-iCiitW^^ki)  [  ’ 

Note  that  /q  is  selected  so  that  (Q  +  WoiA-ih)  <  0.  Substituting  this  result  into  Eq.  24,  the 
new  pseudo-control  becomes 

V2  =  Ides  +  h - =  -  3  im^'ides  (27) 

8  -  (Ci  +  K  2 ki)  S  ~  (C<  +  K  2 ki) 

With  perfect  inversion,  v\  =  Vt  and  v2  =  7.  Thus, 


..RHP  u. 
vO  2K% 


So,  as  desired,  the  system  is  now  decoupled,  however,  the  right-half  plane  zero  is  still  present. 
This  is  the  penalty  for  a  nonminimum  phase  plant.  Figure  1  shows  a  standard  dynamic 
inversion  controller,  while  Figure  2  shows  the  modifications  necessary  to  implement  the 
modified  dynamic  inversion  control  law.  Note  that  both  Figures  1  and  2  are  drawn  assuming 
that  the  D  matrix  in  Eq.  1  is  identically  zero  and  that  the  matrix  CB,  Eq.  9,  is  nonzero  and 
invertible.  A  slight  modification  must  be  made  if  these  conditions  are  not  met. 

Case  2 

Again,  suppose  that  one  transmission  zero,  Q,  is  in  the  right-half  plane  and  would  thus 
result  in  unstable  zero  dynamics.  Now,  let  the  output  zero  direction  be  given  by 


inRHP  q  .RHP 
1  ^0  2 


p 


In  this  case,  multiple  pseudo-inputs  excite  the  bad  zero  state  (the  state  corresponding  to 
the  nonminimum  phase  zero).  In  order  to  stabilize  the  system,  select  one  pseudo- input  to 
stabilize  the  state  associated  with  the  right-half  plane  zero  (£,RHP).  Hence,  let 


Vdesi  +  k{^ 
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Then,  from  the  form  of  wPHP  given  in  Eq.  29,  the  dynamics  of  the  bad  zero  state  become 


cRHP  _  CRHP  ,  „rRHP 
S  —  SiS  +  w0 


Vl 

V2 


—  t P^^P  A-  'inPp^P  i),  -I-  inPP  P  tin  ■  ■  ■  4-  mPPP  11 
~  SiS  +  w0  \vl  +  w0  2^2  +  +  WQ  pVP 


(31) 


=  «BHP  +  •  •  ■  +  (fe  +  k,l;RHP)  +•••  +  «-{ 

Solving  for  £ppp  yields 


ppp 


7?  P 


a..RHP  „f  i  i  „njrLii-r  n.  I  „nrLii-r  •  _i_  _i_  _i_  0, 

W0  1V1  d - r  +  ^0  ;?/desz  +  ^0  i+l^i+l  + - r  ^0 


,PPP 


„MP 


MP 


,RHP 


e 


MP 


^ 

Substituting  Eq.  32  into  Eq.  30  and  simplifying  yields 


(32) 


Vi  = 


(s  -  co  ,  ( wr^vi  +  ■■■+ <“ + <“^1%  +  •  •  • + <“ Vp 


5  -  (Ci  + 


+A)j 


S  -  (C  +  <^*0 


(33) 

Thus,  exact  decoupling  is  not  obtained.  In  order  to  achieve  decoupling,  it  is  necessary  to 
add  another  term  to  cancel  the  effects  of  Vk,k  =  1,  2,  •  •  •  ,p  :  k  ^  p.  Hence,  modify  Eq.  30 
to  the  following: 


vi=^,i+heHF +Y.qj^r 

5S s  +  c> 


(34) 


Substituting  Eq.  34  into  Eq.  31  produces 


CRHP  _  /-  CRHP  I  „rRHP 
S  —  Sis  +  w0 


Vl 

V2 


L  Vp  J 


—  P-pRHP  I  qnPHP  91,  _|_  7/1  RHP  9,  I  I  9 (,iWVP 

—  Sis  +  w0  ]Vi  -t-  w0  2^2  +  r  w0  pVp 


(35) 


= ciy,,p + •  •  ■ + «rp  + ArfHpp + £ 


QijVdesj  )  i  m  m  m  \  „.,RHP 


u#  S  +  Cl  ) 


+  ---  +  wZnr  V. 


pup 


Solving  for  in  Eq.  35  and  substituting  this  result  into  Eq.  34  gives 


V%  Vdesi  hi 


W, 


RHP 


Vdesi  + 

s-((i  +  KHPiki) 


Qij  Vdesj 

s  +  Ci 

3, 3+l 


£ 


(36) 


7  OF  18 


Grouping  like  terms  and  simplifying  Eq.  36  gives 


Vi  = 


(s  -  c i)  Vdesi  +  'Em#  \kiViW0HPj  +  l+( \9ijVdeaj 


s-(Q  +  w^^ki) 


S-  (Ci  +  wffHPiki) 


(37) 


For  complete  decoupling,  the  second  set  of  terms  in  Eq.  37  must  be  identically  zero.  Using 
this  constraint  and  solving  for  v3  gives 


vi  = 


Letting 

produces 


(•®  Ci)  QijVdesj 

Is  +  Q)  kiW™^ 

Qij  =  kiw^HPi 


vj  = 


Ci)  Vdesj 

(s  +  Ci) 


(38) 

(39) 

(40) 


As  before,  with  perfect  inversion,  v1  =  Vt  and  v2  =  7.  Then,  in  this  case,  the  expression 
between  the  desired  and  actual  controlled  variables  becomes 


Vl 

r 

V2 

Vi 

— 

Vp 

- 

—  (s— Ci) 
(s+Ci) 

0 


0 

0 


0 

-(s-Ci) 

(S+Ci) 


0 
0 

0  7-(c 

0  0 


0 

0 


S-Ci 


0 

0 

0 

~(5~Ci) 

(S+Ci) 


1 

Vdesi 

Vdes2 

Vdesi 

- 

Vdesp 

(41) 


So,  as  desired,  the  system  is  now  decoupled,  however,  the  right-half  plane  zero  is  still  present. 
This  is  the  penalty  for  a  nonminimum  phase  plant. 


Results 

For  the  system  described  by  Eqs.  1-  5,  one  zero  is  located  in  the  RHP  (see  Eq.  7).  In 
this  case,  w0  associated  with  the  nonminimum  phase  zero  is 


w0  i  = 


-0.000844  18.866 


For  an  initial  design,  assume  that 


w0i  = 


0  18.866 


0  18.866 


(42) 


(43) 
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so  that  case  1  applies.  Thus,  it  can  be  seen  that  only  the  pseudo-input  associated  with  the 
flight  path  angle,  7,  affects  the  bad  zero  state.  Let 


^-'l  Vldes  dtdes 
=  fade,  +  k(,RHF  =  7j„  +  k 


(45) 

(46) 

(47) 

(48) 


Note  that  a  right  half  plane  zero  has  two  effects  on  the  system,  the  first  a  magnitude  and 
the  second  a  phase.  Unfortunately,  the  phase  effects  cannot  be  altered  with  this  technique, 
however,  the  magnitude  effects  can  be  eliminated  by  choosing  k  so  that  an  all-pass  filter  is 
obtained.  Thus,  the  pole  has  been  placed  at  the  mirror  (about  the  juj  axis)  location  in  the 
left  half  of  the  s-plane.  Now,  the  pseudo-input  relationship  becomes 


1  0  Vtdes 

0  s— 1.949  ■ 

s+1.949  J  L  tdes 


With  perfect  inversion,  the  following  closed-loop  system  is  obtained 


(49) 


(50) 


Figures  3  and  4  show  the  velocity  and  flight  path  angle  responses  along  with  the  ideal 
response  found  using  Eq.  50  and  the  control  deflection  time  histories.  A  small  amount  of 
error  exists  and  this  is  directly  related  to  the  assumption  in  Eq.  43.  It  should  be  pointed 
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out  that  prefilters  were  wrapped  around  the  inversion  controller.  For  the  velocity  channel, 
a  simple  proportional-integral  prefilter  was  used,  while  for  the  flight  path  angle  channel, 
a  proportional-integral-derivative  prefilter  was  used.  The  prefilters  were  the  same  for  the 
simulation  runs  of  cases  1  and  2. 

Now,  relax  the  assumption  in  Eq.  43  so  that  case  2  applies.  Using  Eqs.  38,  39,  and  40 
and  selecting  v2  to  stabilize  the  bad  zero  state  yields 


V2  =  V2des  +  KRHP 


A  des  +  kiRHP  + 


Q21  Vtdes 
~s  +  XRHP 


Selecting  q2 1  =  kwRHP1  gives 


_  .s  -  \RHP  ■ 

1  \  R  T-f  P  *tdes 

S  +  \RHP 

s  -  \RHP 

V2~  s-  ( \RHP  +  wRHP2k)ldes 
Now,  the  pseudo-input  relationship  becomes 


Vl 

— 

'  s-XRHP 

S+\RHP 

0 

0 

s-XRHP 

Vtdes 

V2 

S-(XRHP+WRHP2) 

^des 

With  perfect  inversion,  the  following  closed-loop  system  is  obtained 


Vt  ' 

0  1 

s  \s+XRHP )  U 

Vfdes 

7  _ 

n  1  f  s-XRHP  \ 

s  ys-{XRRp+wRHP2} J 

' ides 

(51) 


(52) 


(53) 


(54) 


Figures  5  and  6  shows  the  velocity  and  flight  path  angle  responses  along  with  the  ideal 
response  found  using  Eq.  54.  Notice  that  the  error  between  the  actual  and  ideal  responses  is 
much  less  than  that  seen  in  Fig.  4.  Again,  this  is  directly  attributable  to  utilizing  the  case 
2  work. 


Conclusions 

In  this  work,  a  dynamic  inversion  type  controller  was  developed  for  an  unstable,  non¬ 
minimum  phase  hypersonic  vehicle  model.  The  technique  used  here  allows  decoupling  of  the 
system,  in  the  same  way  that  standard  dynamic  inversion  allows  decoupling.  The  difference 
is  that  the  nonminimum  phase  zero  cannot  be  cancelled  by  an  unstable  pole.  Hence,  the 
nonminimum  phase  zero  is  retained  in  the  closed  loop  and  a  user  selected  gain  is  used  to 
place  the  left  half  plane  pole. 
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Standard  Dynamic  Inversion 


Figure  1: 


Modified  Dynamic  Inversion 


Figure  2: 


Actual  and  Ideal  Responses 


Figure  3: 
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Time  (sec) 
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Actual  and  Ideal  Responses 


Figure  5: 
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Figure  6: 
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